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Abstract 

The ABC flow is a prototype for fast dynamo action, essential to the origin of magnetic 
field in large astrophysical objects. Probably the most studied configuration is the classical 1 : 
1 : 1 flow. We investigate its dynamo properties varying the magnetic Reynolds number Rm. 
We identify two kinks in the growth rate, which correspond respectively to an eigenvalue 
crossing and to an eigenvalue coalescence. The dominant eigenvalue becomes purely real 
for a finite value of the control parameter. Finally we show that even for Rm = 25000, 
the dominant eigenvalue has not yet reached an asymptotic behaviour. Its still varies very 
significantly with the controlling parameter. Even at these very large values of Rm the fast 
dynamo property of this flow cannot yet be established. 



1 Introduction 

We investigate the kinematic dynamo acti on associated with the well known ABC-flow (which 
stands for Arnold, Beltrami, Childress, see lDombre et al.l [U986I1 ). We focus on the highly sym- 
metric, an d most classica l setup : A : B : C = 1 : 1 : 1 . Its dynamo properties have been assessed 
in 1981 bv lArnord et all 11198 ill and it represents sin ce then the prototype flow for fast dynamo 
action. A "fast dynamo' iChildress and Gilbertl [119950 is a flow which achieves exponential mag- 
netic field amplification over a typical time related to the advective timescale and not the ohmic 
diffusive timescale (in which case it is referred to as a "slow dynamo"). It is known Vishikl [| 1989TI 
that exponential stretching of fluid elements is necessary for fast dynamo action. The existence 
of fast dynamos is essential to account for the presence of magnetic field in astrophysical bodies, 
for which the ohmic diffusive time is often larger than the age of their formation. If self-excited 
dynamo action is to generate their magnetic fields, it is therefore essential that it be achieved 
over an advective time scale. The most clas s ical fl ow to exemplify such "fast dynamo" action is 
indeed the ABC-flow. lArnold and Korkinal Iil983ll first investigated the dynamo property of the 
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ABC-flow, originally introduced to investigate Lag rangian chaos. Many developments foll owed, 
which will be discussed in the course of this articl dGalloway and Frischl [119841 Il986lll987ll . 
Most of the recent developm ents in this field involve non-linear studies wit h forcing be- 



[20051, lArchontis etaD ll2007ll . with a few 



longing to the class of ABC fl ows Courvoisier et al.l - ,. ,_ 

noticeable exception ^Alexakis [ 201 1 jT lGallowayl |20 1 2|] . The asymptotic behaviour of one of 
the most classical example of fast dynamo is however still not understood. This motivates the 
following high-resolution linear study. 



2 Numerical method 



We are concerned with the kinematic dynamo problem, for which a solenoidal magnetic field 
evolution is governed under a prescribed flow by the induction equation 



dB 



_ = Vx (u x B Rm 1 V x B) 

dt v ' 



We consider here the ABC-flow ( Arnoldl [|l965 1. Henon U966ll ). which takes the form 



u = (A sin z + C cos y) e x + (B sin x + A cos z) e y + (C sin y + B cos x) e z 



(1) 



(2) 



and restrict our attention to the case where the mag netic field has the sa me periodicity as the 
flow (i.e. 2^-periodic in all directions of space, see lArchontis et al. 1 fl2003ll for extensions) and 
the weight of the three symmetric Beltrami components are of equal strength {A = B = C = X). 

Let us stress again that we also restrict our attention to the kinematic dyna mo problem, in 
which the flow is analytically prescribed and unaltered by the magnetic field ("see lGalloway and Frisch 
tl987|] for an investigation of the stability of this flow). 

The choice A: B : C = 1:1:1 belongs to the largest symmetry class for this kind of flows, 
and has for this reason been the most intensively studied. Howev er, it yields very small chaotic 
regions and is thus possibly non optimal for dynamo action (see lAlexakisI 11201 ill for a detailed 
study of this point). 

The simulations pre sented in this article were p erformed using a modified version of a code 
originally developed by iGalloway and Frischl lll984ll and which uses a fully spectral method with 
explicit mode coupling. 



The original time-stepping used by IGalloway and Frischl [1198411 relies on a Leapfrog scheme 
stabilised by a Dufort-Frankel discretization of the diffusive term. Introducing X to denote the 
discretized diffusion operator, which is local in Fourier space, and NX, to denote the discretized 
inductive term, non-local as it couples neighbouring modes, this scheme can be expressed as 



B n+i = B „-i + 2dt ( N£ ( B *) + ix(B n+1 + B"- 1 )^ , 



(3) 



using a red-black (or Chloride-Sodium) staggering in time and space, see IGalloway and Frisch 
Il98dl . 
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We have implemented two alternative time stepping schemes, in order to assess the stability 
of the temporal evolution at large values of Rm. We used a Crank-Nicholson Adams-Bashforth 
scheme 

B n+i = B „ + dt f±£(B n+1 + B n ) + | N£(B n ) - 1 NZiB 11 - 1 )] , (4) 

as well as a second order BDF discretization 

^B n+1 = 2 B n - B n 1 + dt (£(B n+1 ) + 2 N£(B n ) - ^(B^ 1 )) . (5) 

These two schemes are unstaggered and involve larger memory requirements, still offering the 
same complexity. All schemes are semi-implicit, but retain an explicit marching for the non-local 
term in order to prevent the resolution of a linear system at each time-step. We verified that the 
results presented in this article are independent of the above choices. 

The computing time obviously varies with the control parameter Rm. If all spatial modes 
are computed up to a truncation N, the computational complexity scales like 0(N A ). Assuming 
the asymptotic scaling of the magnetic field length scale, we get N ~ Rm 1//2 and thus expect a 
complexity growing as (9(Rm 2 ). We have therefore derived a parallel version of the code using 
the MPI library and a spectral domain decomposition strategy to tackle larger values of Rm. This 
yields shorter computing time at large resolution. 

The results presented in this article were obtained with numerical resolutions ranging from 
N = 64 for the smallest values of Rm to N = 1024 for Rm = 25000. In all cases we verified that 
the results reported here were unaltered by doubling the resolution. The simulations presented 
here were performed on up to 512 cores. 

It is worth stressing that the quantity V B is obviously preserved by (Q3, and that this essential 
property is retained by the discrete numerical schemes, and thus the magnetic field remains 
solenoidal throughout the simulations. 

We investigate a linear problem and therefore expect that, independently of the initial con- 
ditions, the long time integration will simply reflect the eigenmode with largest growth rate. In 
practice, we used two different sets of initial conditions, either 

B(t = 0)oc(sinz — cosy) e x + (sinx — cosz) e y + (siny — cosx) e z , (6) 

or a random initial condition with a spectrum converging as k 2 for regularity and projected 
numerically to get a non-divergent field. 



3 Modes crossing 



This study is focused on th e 1:1:1 ABC- fl ow, which has a l arge number of symmetries. These 
have been well documente dArnordl [1198711 . IPodviginal IU999I1 . The symmetry group is generated 
by two independent rotations: the cloverleaf rotation T : x — > y ^ z ^ x and a rotation 
of angle n/2 followed by a translation (Rt in the sequel) around one of the three coordinate 
axes. For example Rt y : x —*■ n/2 + z,y —*■ y — n/2,z —*■ n/2 — x. All the other symmetries 
are obtained by combinations of these two rotations. The resulting group of symmetries of the 



3 




ABC-flow contains 24 elements (including identity). We tested our numerical results against all 
the symmetries of the ABC-flow. 

Figure[T]presents the evolution of the maximum growth rate of the magnetic field as a function 
of Rm. Each point on the figure correspo nds to a three-di mensional simulation. We confirm 
growth rates obtained by earlier studies (see iGallowayl [|2012ll for a recent review), and we extend 
the range of investigation from Rm < 1600 to Rm < 25000. The curve has b een validated 
agains t published growth rates using spec tral methodslGallowav and Frisch 1 1986 1. Lau and Finn 
11199311 as well as a finite volume method iTeyssier et al. [12006 1 for which simulations have been 
performed up to Rm = 2000 (Teyssier & Dormy private comm.). 

In addition to the wider extend of Rm variation, our curve also offers a finer resolution 
than previously obtained graphs. This highlights the presence for two kinks in the curve, la- 
belled Rm i and Rm? on th e figure. The first of these occurs in the stable window reported by 
Galloway and Frisch] I11986ll near Rm = 20 and corresponds to Rm! e [24.05, 24.10] . A mode 



crossin g was previously suggested owing to the changes in the eigenfunction symmetry lGalloway 
1120 121 1 . Here we demonstrate this eigenvalue crossing by following both eigenvalues on each side 
of the crossing. In fact whereas time stepping algorithms usually only provide information on 
the dominant eigenvalue, i.e. the eigenmode with largest growth rate, we use it here to get more 
information. Indeed, transient behaviour starting with well selected initial conditions provide 
information on the behaviour of a given mode, even if it is not the dominant eigenmode (see 
Figure |2]). This transient behaviour allowed us to continue the branches corresponding to each 
eigenvalue outside of the region in which they are dominant eigenvalues (see dotted lines and 
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a. 



b. 



Figure 2: Time evolution of the magnetic energy, which highlights the information provided by 
transient behaviours. Left (a) a simulation for Rm = 16.5 (< Rmi) using as initial condition the 
final solution obtained at Rm > Rmi, right (b) simulation for Rm = 28.4 (> Rmi). 



open symbols on Figure [T). 

In the first window Rm < Rm! as argued by Arnol'd^moTd 1 1987 1. we observe that the dom- 
inant eigenmode has all "even" symmetries of the ABC-flow, i.e. it has every combination of an 
even number of Rt as a symmetry, and is antisymmetric otherwise (the solution changes sign by 
the corresponding transformation). In the second window Rm > Rmi, we observe numerically 



that al l the above symmetries and ant i-symmetries disappear, as pointed by iGalloway and Frisch 
1 1986ll . However, Jones et al. [ 2012 1 are currently using a decomposition of this mode in three 



components each satisfying different symmetries. 

As noted by earlier authors, the dominant eigenvalues are complex, leading to oscillations of 
the energy, visible on Figure |2] (particularly on the first part of Figure |2k, as the period of oscil- 
lations is elsewhere very short compared to the time extend of the plot). The imaginary part of 
the dominant eigenvalue co = Im (A) can thus be directly determined from these time series. The 
graph w(Rm) is displayed on Figure|3] As explained above, not only do we display the dominant 
eigenvalue (solid line and symbols) but we are also able to follow each mode past their region 
of selection (dotted lines and open symbols). One can note that, as Rm increases past Rmi, the 
imaginary part of the dominant ei genvalue jumps discontinuo usly from co ^ 0.53, corresponding 
to the first windo w identified by Arnold and Korkinal IU983I1 . to w ^ 0.13 corresponding to the 
second window of IGalloway and Frisch] [U986I1 . This discontinuous jump in the pulsation high- 
lights the eigenvalue crossing occurring at Rm^ There again, transient behaviours were used to 
obtain the open symbols. 



4 Oscillatory dynamics 



In order to improve our understanding of the oscillatory dynamics, we introduce a phase space 
for this linear system. We rel y for Rm < Rm i on two vectors, corresponding respectively to the 
dominant Ikl = 1 contributio nArnordl [|1987] 



bi = (sin(z) — cos(y)) e x + (sin(x) — cos(z)) e y + (sin(y) — cos(x)) e z , 



(7) 
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Figure 3: Plot of the imaginary part of the eigenvalue as a function of Rm (using logarithmic 
scale in the x-axis). 



and to Ik I 



b2 = (sin(y) cos(z)) e x + (sin(z) cos(x)) e y + (sin(x) cos(y)) e z . 



(8) 



These two vector fields satisfy all symmetries of the realised eigenmode for this parameter regime 
(Rm < Rmj). 

Simulations in this regime rapidly reach an asymptotic behaviour starting with bi as an initial 
condition. 

We construct the phases by introducing 



X = <B • bi> exp - ^ , and Y = <B b 2 > exp~ 



(9) 



where the growth rate cr is a function of Rm (see figured])- As the governing equations are linear, 
the exponential damping is here essential in order to introduce a limit behaviour. The quantities 
X and Y are presented on figure for Rm = 1 1 .5. The undamped trajectory of the system is also 
represented using dashed lines, and directly illustrates the exponential growth of the dominant 
mode. The damped trajectory evolves toward the equivalent of a stable limit cycle. 

The exponential damping on X and Y provides us with the equivalent of a non-linear dy- 
namical system. This explains why we report below phenomena that are usually associated to 
non-linear dynamics. 

In the second window, Rm > Rmi, the two modes we selected (corresponding to the lowest 
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Figure 4: Phase space diagrams for Rm < Rrri! (a) and Rm > Rrri! (b). 



k component of the realised mode) are 

h\ = (— s'm(y) — cos(z)) e x + (sin(z) + cos(x)) e y + (— sin(x) + cos(y)) e z 
b' 2 = (sin(y) — cos(z)) e x + (sin(z) — cos(x)) e y + (sin(x) — cos(y)) e z . 



(10) 
(11) 



These correspond to two components of the general family reported earlier for this mod dGalloway and Frisch 
rtl986ll . These two modes involve symmetries (though none in common), but these are not rele- 
vant here since these are not verified by the full eigenmode. 

The resulting orbits of X' = (B ■ b'i) exp _lTf and Y' = (B • b' 2 ) exp~°"' are represented on 
figure 0b for various values of the controlling parameter Rm > Rmj. The represented quantities 
are not arbitrarily rescaled. Only the exponential damping has been applied, and all cases have 
been started with the same initial condition (involving random, but divergence free, fluctuations). 
The oscillating nature of the dynamo for Rm < Rm 2 is clearly illustrated by the limit cycle. For 
Rm > Rm 2 the oscillations disappear and the dynamo mode therefore becomes a fixed point in 
the [X', Y'] plan. For clarity, we suppressed the trajectories that lead to the limit cycles or the 
steady solutions on figure 



5 Eigenvalues coalescence 

As the magnetic Reynolds number is further increased, a second kink in the growth rate is ob- 
served on Figure [TJfor Rni2 e [215.0, 215.4] . This second accident, however, does not correspond 
to a change of dominant eigenvalue, but instead to an eigenvalues coalescence. The strategy high- 
lighted above to follow secondary modes is inefficient here, indicating that there is no significant 
change in the dominant eigenmode. 

Figure |3] reveals that the behaviour of the imaginary part of the eigenvalue is very different 
near the second kink. Instead of the abrupt jump reported at Rm 1; the pulsation continuously 
(but not smoothly) tends to zero as Rm approaches Rm 2 and vanishes for Rm > Rm 2 . 

The lack of oscillati o ns at large Rm is a well known characteristic, it was already noticed 
by iGalloway and Frisch [Il984ll for (Rm > 400), although they could not assess whether the 



period of oscillations was simply increasing with Rm or the eigenvalue had become purely real. 



Lau and Finnl (1199311 suggested that this could be associated with a mode crossing, a new mode 
with purely real eigenvalue taking over above Rm 2 . 

We show here that the imaginary part of the eigenvalue indeed vanishes for Rm > Rm 2 , and 
that this corresponds to the coalescence of two complex conjugate eigenvalues on the real axis. 
The coalescence yields the kink in the evolution of the real part of the eigenvalue. 

The simplest mathematical model for a complex conjugate eigenvalue coalescence on the 
real axis corresponds to a situation of the form 



A+ = or(Rm) + y/3(Rm) , (12) 

where a and are differentiable real functions of Rm. A negative /? (for Rm < Rm 2 ) yields two 
complex conjugate modes, and thus oscillations of the magnetic energy. As (5 becomes positive 
(for Rm > Rm 2 ), the eigenvalues are purely real and the /? term now contributes to the real part 
of the eigenvalue A + offering the largest growth rate. 

Figure |5^ presents a detailed view on the variation of cr = Re (A) and co = Im (A) close to 
Rm 2 . Defining cr 2 = cr(Rm 2 ) we plot cr(Rm) — <r 2 and -w(Rm). It is clear that the kink in cr is 
concomitant of the vanishing of co. 

Let us now form on Figure |5t> the quantity F = <x(Rm) — cr 2 — w(Rm). The square-root be- 
haviour of F near Rm 2 is obvious. Assuming that the above model (PT2l) is correct, F corresponds 
to a — a + sign (J3) \f\j3\, where ao = a(Rm 2 ). We can note on Figured for Rm < Rm 2 that 
a — ao remains small compared to variations in f3. The quantity sign (Rm — Rm 2 ) F 2 therefore 
offers a good approximation to J3 and should be differentiable at Rm 2 . 

More formally, assuming that a and ft are regular functions of Rm, we can write a finite 
expansion of the form 



a = 



= a + ori(Rm - Rm 2 ) + ar 2 (Rm - Rm 2 ) 2 + • • • , (13) 
J3 = ft(Rm - Rm 2 ) + /? 2 (Rm - Rm 2 ) 2 + • • • , (14) 

with fii > 0. The quantity sign (Rm — Rm 2 ) F 2 can be written at the lowest orders in Rm — Rm 2 

sign (Rm — Rm 2 ) F 2 = sign (Rm — Rm 2 )ar 2 (Rm — Rm 2 ) 2 

+ sign(Rm-Rm 2 )a 1 v ^ 1 (|Rm-Rm 2 |) 3/2 (15) 
+ /? 1 (Rm-Rm 2 ). 

This development implies that sign (Rm — Rm 2 ) F 2 is differentiable at Rm 2 . Figure |5t> clearly 
illustrate this property on the direct numerical simulation. 

Another insight on the nature of this transition can be gained from the "phase space" intro- 
duced in figure m>. As noted above a remarkable feature is that the steady solutions obtained 
after the coalescence lies on the ellipse described by the limit cycle shortly below the coales- 
cence (on the figure Rm = 220). This behaviour which is similar to that of an excitable system 
(such as a pendulum s ubject to a constant t orque), has recently been observed in experimental 



dynamos measurement aRavelet et al.l [1200811 . The phase space was then constructed using two 



components of the magnetic field at a given location (a probe). 
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Figure 5: Eigenvalues coalescence for Rm close to Rm 2 ^ 215. Both the real part (written as 
cr(Rm) — cr 2 ) and the opposite of the imaginary part (i.e. -w(Rm)) are represented (a). The sum 
F of both (b), and sign (Rm — Rm 2 )F 2 illustrate the continuity and regularity in the functional 
form of the eigenvalue (see text). 



Such behaviour is reminiscent of a saddle-node bifurcation. In these systems, the dynamics 
becomes increasingly slow on the cycle as the system approaches the state at which the sad- 
dle and the node will collapse. In order to assess this property in our system, albeit linear, 
we introduce the angle 9 of the system over a unit circle described through the orbit, so that 
(cos(0), sin(0)) = (X', Y')/ \/X' 2 + Y' 2 . The time evolution of 9 with increasing values of Rm is 
presented in figure |6] The system clearly spends an increasing amount of time as Rm approaches 
Rm 2 near the angle at which the stable solution will occur for Rm = Rm? . Such behaviour coul d 



be des cribed by a simple phase dynamics, e.g. IPetrelis et al.l [|2009| 1. iGuckenheimer and Holmes 
lll997ll . 

It is interesting that these approaches of non-linear dynamics can cast some light on the be- 
haviour of kinematic dynamos. A similar occurrence of a saddle-node transition for a kinematic 
dynamo numeric al model of the above mentioned VKS experiment has indeed been reported in 



Gis singer [2009], with an expression of the form (fT2l . 



The fact that the fixed points in the (X', Y') plan describe the previously existing limit cycle 
is a strong indication that there is at first no significant change occurring in the structure of the 



dominant eig enmode after the eigenvalues coalescence. Indeed the "double cigars" structure (see 
h| ll2000ll '). associated to the oscillations for Rm e [Rm 1? Rm 2 ] is preserved once the growth 



Dorcl 



rate has become steady, Rm > Rm 2 (see Figure [7]). 



6 Asymptotic behaviour 



We have finally increased the control parameter in the range 2 15-25000. Despite the fact that the 
largest magnetic Reynolds number tackled in this study is roughly 15 times larger than earlier 
results, the growth rate has not reached an asymptotic value yet. The growth rate obtained for 
our largest Rm is very close to 0.1 and appears to be stil l significantly i ncreasing with Rm. 

The 1:1:1 ABC-flow has also been considered by IGilbertl II 199211 using maps in a limit in 
which the diffusivity is formally set to zero. This approach has yield growth rate of 0.04 - 0.05, 



9 



71 



7t/2 





Rm = 200 




Rm = 210 




Rm = 215 




Rm = 220 



-71/2 



1000 



2000 



Figure 6: Evolution in time of the projection of the dominant mode on the two large scale com- 
ponents, represented via the angle 6, such that (cos(6>), sin(0)) = (X', Y')/ \/X' 2 + Y' 2 . 




Figure 7: Eigenmode obtained for Rm = 434.2 (> Rm 2 ). An isosurface of the magnetic energy 
is represented. To highlight the symmetry, it is coloured according to the value of B x , from blue 
(negative) to red (positive). The corresponding eigenvalue is purely real, yet the double cigars 
structure of the field remains clearly visible. Magnetic field lines are also depicted on the right 
plot. 
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so much smaller than the value achieved by our direct numerical simulations at Rm = 25000. It 
is therefore not unplausible to anticipate that the behaviour of cr(Rm) above 25000 will not be 
monotonic and cr will probably decrease again. 

Another indication i s provided by the la rgest Lyapunov exponent of the flow, which is ap- 
proximately 0.055 (see Galanti et al.1 ||l992|l ). Owing to the lack of regularity of the field in 
the limit of large Rm numbers, the larges t Lyapunov exponent howev er does not provide an 
upper bound on the asymptotic growth rat aChildress and Gilbert [19951. An upp er bound can 
be sought by considering the topological entropy h top (see Finn and Ott 1 1988a For steady 
three-dimensional flows the topological entropy is equal to the line stretching exponent hn ne (see 
Childress and Gilbert! [ Il995h . which can be estimated for the 1:1:1 ABC-flow to be h\ ine ~ 0.09. 



This provides yet another indication that the curve cr(Rm) must decrease for larger val ues of Rm. 



A plausible scenario, suggested by the behaviour of submodes as investigated by lJones et al. 



II2OI2I1 . is that two complex conjugate eigenvalu es may emerge agai n at larger Rm. This is 
often observed in saddle-node bifurcations (e.g. IRavelet et al.1 [2008|]). It would result in the 
reappearance of the oscillations, and an abrupt decrease of the growth rate (the counterpart of the 
increase observed at Rm 2 ). This would deserve further study. 

The asymptotic behaviour of the 1:1:1 ABC-flow is thus not yet established. It is at the mo- 
ment, despite the high resolution simulations presented here, impossible to assess its asymptotic 
growth rate. It is not even possible to rule out the possibility of an eventual decay of the growth 
rate to zero at very large Rm. 



7 Conclusion 



We have investigated using high resolution direct numerical simulations the behaviour of the 
1:1:1 ABC-dynamo. We have shown that the two dynamo windows identified for this dynamo 
are associated with a change of dominant eigenvalue. We have identified a second kink in the 
growth rate as a function of Rm and shown that it corresponds to an eigenvalue coalescence and 
the end of the oscillatory nature of the solutions. Finally, even at very large values of Rm, we 
show that the growth rate is still strongly varying and not monotonic yet. 

Relaxing the requirement of a fully three-dimensional flo w and allowing for time depe ndence, 
other models for fast dynamo actions have been obtained by iGalloway and Proctor! [ |1992|l . with a 
velocity depending only on two coordinates. The time dependence ensures exponential stretching 
at least in this plane. The induction equation is then separable in the z direction, allowing faster 
numerical integrations. The asymptotic limit of large Rm appears easier to reach for such flows. 
Other studies involved time dependence of the flow Otanil [U993I1 and some hint at a possible 
resonance phenomenon Dormy and Gerard-Varet|[t2008ll . 

Finding a good example of fully three-dimensional flow that acts as a fast dynamo remains a 
challenging problem. The most classically given example remains the 1:1:1 ABC-flow. This 
unexpectedly rich behaviour of the 1:1:1 ABC-dynamo at very large Rm, highlighted in our 
study, deserves further in vestigations. It is most likely associated with the fact this flow yields 
very small chaotic regions lGallowayl B2012fl . 
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